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We present a theoretical framework to describe the collective emission of light by entangled atomic 
states. Our theory applies to the low excitation regime, where most of the atoms are initially in the 
ground state, and relies on a bosonic description of the atomic excitations. In this way, the problem 
of light emission by an ensemble of atoms can be solved exactly, including dipole-dipole interactions 
and multiple light scattering. Explicit expressions for the emitted photonic states are obtained in 
several situations, such as those of atoms in regular lattices and atomic vapors. We determine the 
directionality of the photonic beam, the purity of the photonic state, and the renormalization of the 
emission rates. We also show how to observe collective phenomena with ultracold atoms in optical 
lattices, and how to use these ideas to generate photonic states that are useful in the context of 
quantum information. 



I. INTRODUCTION 

In the last years, the field of atomic, molecular and op- 
tical physics has witnessed an impressive advance in the 
development of setups to trap atoms under different con- 
ditions, like for example, ions in electromagnetic traps 
and ultracold atoms in optical lattices. Furthermore, 
the quantum state of these systems may be engineered 
by performing quantum operations such as as quantum 
gates between ions [l| , or the excitation of neutral atoms 
under the dipole blockade @, Q . In this way one can cre- 
ate deterministically collective entangled states, like the 
completely symmetric states with a single excited atom 
(|W)-states) Q. In addition, those are systems where 
atoms can be coupled to light in a very controlled way. 
Since some of the atomic entangled states which may 
be created in those setups play an important role in the 
description of the interaction of atomic ensembles with 
light, an important question arises, namely, can we use 
our control on the states of trapped atoms to generate 
useful quantum states of light? If so, which are the prop- 
erties of such states in terms of photon directionallity, 
purity, or photon entanglement? 

A preliminary study of these ideas was recently pre- 
sented in [j| . The understanding of this problem requires 
a theoretical framework to describe the emission of light 
by collective atomic states under a variety of trapping 
conditions. Related problems have been indeed a subject 
of investigation since the seminal work by Dicke [6| . Nev- 
ertheless, the experimental systems that we have in mind 
share a few peculiarities which demand a new approach 
for their description. First, ultracold atoms are trapped 
in ordered arrangements such as Coulomb crystals or op- 
tical lattices, where the distance between atoms, do, is 
comparable to the wavelength of the light, A. This situa- 
tion is very different from Dicke superradiance, where A 
is larger than the size of the whole system. It also differs 
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from the case of crystals, where A 3> do as considered, 
for example, in 7]. Also, we study a situation in which 
atoms are initialized in a collective state, which is not 
necesarilly created by the absorption of a photon. Thus, 
rather than the more traditional description in terms of 
light scattering, we need a theoretical model of the map- 
ping between atomic collective excitations and photons. 

In its more general form the above described situa- 
tion poses a very complicated many-body problem. A 
crucial simplification is achieved by considering the low 
excitation sector of the atomic Hilbert space, which in 
the Holstein-Primakoff (HP) approximation can be de- 
scribed in terms of bosonic spin- waves [|| . This approach 
has been used in previous works, for example, to s tudy 
Dicke superradiance @, slow propagation of light [Tof . 
and atom- light interfaces with atomic vapors Also, 
the emission of light by ensembles of harmonic oscilla- 
tors was studied in [12j | , although in a different regime of 
trapping conditions than those considered here. 

In this paper we make use of the HP approximation 
to describe the collective emission of light as a mapping 
between spin-wave excitations and photons. For atoms 
placed at fixed positions this mapping is a Gaussian com- 
pletely positive map [l3|, and we present a method to 
get explicit expressions of the photonic modes into which 
light is emitted. Furthermore, we extend this formalism 
to study the case of atomic vapors. We obtain the fol- 
lowing results: (i) For atoms trapped in a regular lattice, 
there is a regime in which photons may be emitted in a 
collimatcd beam, which requires that A > 2do- In this 
regime, there is a renormalization of the emission rates, 
leading to a classification of the low-excitation atomic 
Hilbert space in terms of superradiant and subradiant 
spin-waves. Superradiant states decay with a rate that 
is enhanced by a factor x, which depends on the dimen- 
sionality of the lattice. In ID and 3D we determine the 
values xid oc A/do and \3D oc (A / do) 2 (L / do) , respec- 
tively, with L the length of the lattice. This effect is 
related but not equivalent to Dicke superradiance. (ii) In 
the case of atomic vapors the collective characeter in the 
emission of light is determined by Xon oc N(X/L) 2 , where 
N is the total number of atoms. The directional regime 
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requires x C n 2> 1, and the emission rate is enhanced by 
Xcn- (hi) Some of these effects, like the renormalization 
of the emission rates and the directionality, could be ob- 
served in a relatively simple experiment with atoms in 
optical lattices, (iv) By making use of coherent effects, 
photonic entangled states could be generated by trapped 
ions or neutral atoms. Photons that are generated in this 
way could be collimated and in a pure state, and thus be 
useful in the context of quantum information. 

II. THEORETICAL FRAMEWORK 

Our first task is to describe the collective emission of 
light by an ensemble of atoms. First, we focus on the 
situation in which atoms are placed at fixed positions. 
At the end of the section we discuss the effect of atomic 
motion. 

A. The atom-photon map in the 
Holstein-Primakoff approximation 

We consider an ensemble of -/V atoms trapped by har- 
monic potentials with trapping frequency v, and inter- 
nal levels forming a A-scheme, see Fig. [TJ Two ground 
states \g), \s), are coupled by means of a laser with Rabi 
frequency S1l and wavevector kL, through and auxil- 
iary level |a). We define atomic operators Oj = |g)j(s|, 
cr^ 5 — \s)j(s\, which fulfill the commutation relations, 

([a 3l a+])=5 u {l-2(af)). (1) 

Under the condition that the excitation probability of 
each atom is low, (<t^ s ) <C 1, we replace atomic operators 
by HP bosons, <7j — > bj. 

The HP approximation allows us to recast the atom- 
light Hamiltonian as a bosonic quadratic Hamiltonian. 
After the adiabatic elimination of the upper level |a), 
our system is described by 

H = i?o + Him + H mot . (2) 

Hq = ^ k ajfea^ak, with uok — ck. We consider, for the 
shake of clarity, a scalar model for the electromagnetic 
field, since our conclusions do not change when includ- 
ing the dipole pattern, as we show later. However, the 
photon polarization can be included straightforwardly in 
our formalism. H\ m is the atom-light interaction Hamil- 
tonian in the rotating wave approximation [14j | (we set 
71=1), 

#im = ^.g fc fe+a k e' i(k - kL ^ +lWLt +h.c. , 

9k = 2A2^ a - (3) 

V is the quantization volume of the electromanetic field, 
eo the vacuum permitivity, and d ga the dipole matrix el- 
ement of the |ff)-|a) transition. H mot is the Hamiltonian 



describing the motion of the atoms during the emission 
process, whose energy scale is given typically by the trap- 
ping frequency v. Eq. ((3]) is obtained under the assump- 
tion A 3> n^,h i n, with n the mean vibrational occupa- 
tion number. In the resonant case, A = 0, the validity 
of the adiabatic elimination of \a) requires that the de- 
cay rate from \a) to \g), r ag , satisfies r ag 3> f^L, vn^ and 
a similar coupling as ([3]) is obtained. The spontaneous 
decay of \a) back to |s) is neglected in the forecoming 
analysis. This is justified either when collective effects 
enhance the \e)-\g) channel, or by choosing \e)-\g) to be 
a cycling transition, see Appendix |Dj We deal first with 
the simplest case, in which atoms are placed at fixed po- 
sitions and the motion of the atoms can be neglected. 
This is a good approximation in the Lamb-Dicke regime, 
(k^xo) 2 (2n + 1) <C 1, where xq is the ground state size in 
the harmonic trap. The effect of the motion is considered 
in the last subsection. This work relies on the observation 




lg> l*> 

FIG. 1: Lambda scheme. 

that Ho + H\ m generates a beam-splitter transformation 
between atomic states and photons. The evolution of the 
atom-photon system can thus be understood in terms of 
a Gaussian completely positive map [l3| . In the low exci- 
tation regime, any atomic state can be expressed in terms 
of bosonic excitations, 

|*>at= £ *„ N ^...^^|0) at . (4) 
ni,...,n^v 

We consider the following initial state, 

|^(0)) = |*) at |0) phj (5) 

where |0) p h is the photon vacuum. Our goal is to find 
the photon state \$) p h at a time t longer than the atomic 
decay time, 

\m) =w(t)iv(o)) = io) at |$) ph , (6) 

where U(t) = e~ tHt , and H is the total atom-system 
Hamiltonian. Together with the beam-splitter form of 
© this implies that the problem is reduced to finding 
the exact form of the transformation, 

U(t)btU(t? = $> ik (t)a+ + (7) 
k ; 
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Taking into account that U(t)^\v&c) — 0, with |vac) = 
|0)at|0) p h, then gjk(t) can be found by using the relation, 

9jk{t) = (vac|a k W(t)6+|vac) = (vac|a k (i)6+ |vac), (8) 

where a^(t) is the photon operator in the Heisenberg pic- 
ture. Also, 

hji{t) = (vac\btU(t)bj\vac) = (vac|6 ; (i)6+|vac), (9) 

which shows that hji (t) = for t longer than the atomic 
decay time. gjk(t) is thus the only interesting term, since 
it determines the photon mode into which light is radi- 
ated. Its determination can be readly done by means 
of the Heisenberg equation of motion for the photonic 
operators, which yields 



(10) 



^ Te -i(k~k L )r 1 +i(wfc-WL)r 



(vac|6/(r)6~'"|vac). 

This equation is our starting point for the exact deter- 
mination of the photonic modes. 



is the single atom decay rate. The inclussion of the pho- 
ton polarization would change the spatial dependence of 
the couplings Jy. However, as we show later, the collec- 
tive phenomena would be the same. Note that we are not 
including the single atom Lamb shift, which may be sim- 
ply absorved into the laser frequency wl- Since Eq. (jllj) 
is quadratic in bosonic operators, it is readly solved by 
defining eigenmodes which diagonalize the atomic quan- 
tum dynamics in the low-excitation limit, 



bj = 2 ,-Mjnbn, 

n 

(M- l JM) nm = J n S n , m . 



(15) 



The matrix is not hermitean, and thus canonical com- 
mutation relations are not conserved, [b n , b+] ^ S n . m . 
However, the evolution of averages takes a simple form 
given by 



(h(r)) =^M ln e- J » T (6 n (0)). 



(16) 



Thus, the spin-wave dynamics is governed by the eigen- 
values J n . The latter contain the collective decay rates, 
r„, and the collective energy shifts, 



B. Photonic modes 

To find the mapping between atomic and photonic 
modes, we need to solve the master equation, to calculate 
the time evolution of the atomic correlator in Eq. ((TO]). 
In general, this is a difficult many-body problem, but in 
the bosonic limit it can be described exactly. In terms 
of HP bosons, the master equation which describes the 
atomic dynamics reads 15[ , 



dp 
~dt 



= ^2 (Jijbjpbf - J ij bfb jP + h.c.) , (11) 
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where the Jjj include multiple light scattering and dipole- 
dipole interactions, 

/>oo 

./,., V.^ / e *»*-<K)r+iO'-**)l'*-'j)dT. (12) 

k Jo 

This expression is evaluated by using the identity 
J dTe %u)T = ttS(u>) + iV(1/oj), which yields the result, 



Ji i 



p e -ik L (ri-rj 



sin(A: L |r f — ry|) . cos(fc L |r* - Vj \ ) 



ki,\r, 



where 



3tt \2AJ e c 3 ga 



(13) 



(Mi), 



(14) 



(17) 



r n =2Re(J n ), A„ = 2Im(J„ 



Conservation of the trace under the transformation (jT5 
leads to the following sum rules, 



(18) 



Note that whenever collective effects induce a renormal- 
ization of T n , the sum rule implies the existence of super- 
and subradiant states. 

By application of (|15U16[) . and the quantum regression 
theorem, we determine the two-time atomic average in 
(TIT)]) . In the limit that t > 1/T n we get 



9 At) = ze-^gk 



. (19) 



which yields the explicit form of the atom-photon map- 
ping. By means of this relation it is possible to de- 
termine the many-photon state emitted by any initial 
atomic state like (|4]). 

To get insight of the characteristics of the photonic 
states emitted by the collective atomic states, we focus 
from now on, on the mapping to a single photon. To 
clarify the notation, let us define |^ n ) a t 5 the n-spin-wave 
state with a single excitation, 



(20) 
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M has to be included due to the non-hermiticity of Jy. 
We determine |$ n ) p h, the single photon state into which 
|*„) at is mapped, 

|$n(*))ph = J2^ n ^ a i^P h ' 



i — > — s — • (21) 



*Y i(u k - u L ) - J n ' 



Although these results allow one to solve exactly the 
problem of collective emission of light including the ef- 
fects of reabsorption, aditonal insight can be gained by 
considering the case of a system with periodic boundary 
conditions. This will be a good approximation for a finite 
system, provided that the number of atoms in the volume 
is much larger than in the surface, that is, L/do 3> 1. In 
this case we get the matrices 



M 3n - e 



The spin-wave state n is then defined like 
|* n ) at = -L^e^6t|0) a 



(22) 



(23) 



The vector K„ is the momentum of the collective atomic 
state. The resulting photonic mode is defined by the 
following expression, 



0n,k(£) = igke 



-i(k-k L -K n )n 



i(u k - Wl) - Jn 



(24) 



such that collective effects and dipole-dipole interactions 
enter through the dependence of the collective emission 
rate J n on the mode number, n. Note that by using 
periodic boundary conditions, the matrix M. define an 
unitary transformation, and thus the sets |^ r n ) n=1 N , 
and |$ n ) n _j jv) form an orthogonal basis of spin-waves, 
and photonic modes, respectively. 



C. Angular photon number distribution 

We determine now the properties of the emitted pho- 
tonic modes. In particular, let us define 



v r°° 



(25) 



the average photon number per solid angle. Consider an 
initial atomic state with a single excitation, 

Mat = 5>i#|0> rt . (26) 

3 

The emitted photon distribution is 

3,j' 

hi' (") = I 9-*9fMk 2 dk. (27) 



Upon substitution of (| 19|) in the expression for Ijji(Cl), 
and under the condition that Y n /c <ti 1/L (see Appendix 
|A|. we get, 



-k L )r, 



(28) 



where uq is a unit vector pointing in the direction of the 
solid angle Q. In spherical coordinates, Q is determined 
by (0,4>), such that uq = (sin 8 cos 4>, sin 9 sin (j>, cos 9). 

The general recipe for calculating the photon distribu- 
tion probability involves the following steps: (i) Calcu- 
late the coefficients of the master equation and find the 
eigenvalues and eigenvectors of Jy. (ii) Use Eq. (|28l) to 
calculate ly. (hi) Use the latter to calculate the emis- 
sion spectrum with the wavefucntion of any given initial 
atomic state expressed in terms of bosonic spin- waves. In 
the planewave approximation (|22p we can get close ex- 
pressions for the photon distribution. We focus again in 
the single photon case, and define / n (f2) as the photon 
number distribution corresponding to the photonic mode 
emitted by |4 , n )at- By using Eq. (|28|) we get 



(k n -k L -K„)(r,-r,,) 



4tt r„ 



(29) 



3,3' 



This result has a clear interpretation in terms of inter- 
ference of light emitted by the atomic system [l6j]. Eq. 
(|29f not only allows one to calculate the angular emission 
probabilty, but also, due to the normalization condition 

/ j n (n)dn = i, 



h-= f dn— W(kn-k L -K„)(r 3 - 

r / 4?r ^ 

3,3' 



(30) 



This expression provides us with a simple way to deter- 
mine the collective rates under the planewave approxi- 
mation. 



D. Effects of the atomic motion 

Finally, we discuss the effect of the atomic motion on 
the light emission. In the most general case, the inclusion 
of the motional degrees of freedom poses a very compli- 
cated problem which goes beyond the scope of this work. 
Two time scales determine this problem. First, T mot , the 
time scale of the motion of atoms in the trap. In the case 
of trapped particles, r mo t = with v the trapping 

frequency. We can extend this discussion to the case of 
atomic vapors, and consider that in this case, T mot = L/v, 
with L the length of the sample, and v the atom veloc- 
ity T mot is to be compared with the radiative decay time, 
1/f , or more specifically, the set of collective decay times 
1 /r„. Based on the comparison between these time scales 
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we define two limits in which the application of our the- 
oretical framework is particularly straightforward. 

(i) Slow motion limit, r mot 3> l/r n . Since the emission 
process is much faster than the motion of the particles, 
we can assume that atomic positions are frozen. The 
system is in the initial state 



|V(0)) = |tt)at|0)ph|*m)mot, 



(31) 



where x P m (ri, . . . , r/v) is the initial wavefunction in terms 
of the atomic positions, which in this limit does not 
evolve during the emission time. The atom-photon map- 
ping can be still applied to this system by solving it for 
each value of the atomic positions, 



|*>at|0)ph- |0)at|$ r 



/ph, 



(32) 



where | < E ,ri, '"' rw ) p h is the photonic state obtained un- 
der the assumptions that ions are located at positions 
ri, . . . rjv- Any photonic observable, O, is then obtained 
upon averaging with respect to the wavefunction \& m , for 
example, 



(O) 



E ■ 

ri,...,rjv 



4> r 



|OI$ r 



|*(n 



,r N ) 



(33) 

This method can be readly extended to the case of a 
mixed motional state. 

(ii) Fast motion limit, r mot <C l/r n . In this case, we 
are assuming that trapped atoms move along the sample 
in a time that is smaller than the emission time. This 
case is particularly relevant, since it describes hot atomic 
vapors. We notice first that condition r mot -C 1/E„, 
implies that the atomic positions are not correlated with 
the atomic operators in Eq. ([ID)). Thus, we can describe 
the atomic radiative decay indepently of the evolution of 
atomic positions r; in (|10p . This can be done by using a 
master equation with averaged coefficients, 



(k-k L )( ri -r 3 ) 



(34) 



In the case of an harmonic trap, are position opera- 
tors, and (. . . ) mo t is the average with the atomic motional 
state. This situation can be extended to describe hot 
atomic vapors, by replacing the atomic positions, r,, by 
a set of random variables with a given probability distri- 
bution /O m (i"i, . . . , r/v), and performing the corresponding 
average to get J^-. 

Once the atomic dynamics is solved, our results on the 
atom-photon mapping can be used to describe the emit- 
ted photons. A general atomic state is mapped now into 
a mixed photonic state. For example, consider the atomic 
state |\E' m )at, defined by Eq. ([20f . which is obtained by 
the diagonalization of the master equation with (|6ip . The 
atom-photon mapping yields the following photon den- 
sity matrix after the emission process, 

P Ph (t) = X)<^(*Xk'(*)>nurtOt|0)p h <0|o V , (35) 
k,k' 



and this expression is easily generalized to the multipho- 
ton case. Note that the solution of the fast motion limit 
seems similar to the case of slow motion. However, the 
crucial difference is that in the fast case, we are allowed 
to solve the radiative emission problem, and to perform 
subsequently the spatial average in ([33)1 . We will study 
in more detail this situation later in the case of an the 
collective emission properties of atomic ensambles. 



III. ATOM-PHOTON MAPPING IN A SQUARE 
LATTICE 



The situation in which atoms are arranged in a crystal 
is found in experimental setups such as ultracold atoms 
in optical lattices and Coulomb crystals of trapped ions. 
The results presented in the previous section are applied 
here to study the collective emission process in these sys- 
tems. We obtain analytical results by using the plane- 
wave approximation. 



A. General Discussion 

Let us study for concreteness the case of one (ID) 
or three dimensional (3D) square lattices, although our 
results are easily generalized to different lattice geome- 
tries. Assuming periodic boundary conditions, the al- 
lowed wavevectors are 



K, 



2tt \ - n c 



-a. 



(36) 



In the ID case we consider that the atom chain is aligned 
in the z direction. Thus, in (I36|) and the forecoming ex- 
pressions, a runs over a = z in ID, and a = x, y, z in 3D. 
a is a unit vector in the direction a, and N a is the num- 
ber of atoms along a. Considering, for concreteness, the 
case of even N a , each wavevector is determined by the 
set of integers n a = —N a /2, . . . , N a /2 — 1. By applying 
Eq. ()29p we determine the photon number angular dis- 
tribution for the photonic state emitted by a spin-wave 
excitation with momentum K n , 



7 n (Q) 



i r 



(37) 



sin 2 ((fc L uB-kg-Kg)d JV a /2) 
sin 2 ((/c L u«-k£-Kg)d /2) 



In(Q) shows a series of diffraction maxima at solid an- 
gles f2 at which the sin 2 -function in the denominator van- 
ishes. Note that by including the photon polarization, we 
would have got an additional function of O multiplying 
the photon distribution I n (Q), which would correspond 
to the single atom dipole pattern. The latter would in- 
duce the suppression of diffraction peaks, if they are in 
a direction forbidden by the dipole pattern. Since we 
are specifically interested on collective effects we do not 
consider this effect in the discussion that follows. 
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To get a quantitative description of the photon distri- 
bution wc notice first that 



1 sm 2 (Nx) 
N sm 2 (x) 



J2fN{x-mir) if (JV»1), (38) 



where the function fx (x) describes the shape of each of 
the diffraction peaks, 

f N (x) = Nsinc 2 (Nx), -tt/2 < x < tt/2, 

/at (x) = 0, otherwise. (39) 

Thus in the limit N a 3> 1, one can approximate the 
emission probability as a sum over Bragg scattering con- 
tributions, 



7 n (0) = 



4 m] («) 



i r 



(40) 



fs. [— r ~ W J + m Q 7r 



Each term in the sum is labeled by the vector m, and cor- 
responds to a different diffraction peak. The probability 
that the spin- wave n emits a photon in the m diffraction 
peak is given by 



(41) 



Eq. (|40|) has a clear interpretation in terms of momen- 
tum conservation. 7 n (f2) has a maximum whenever there 
is a value of m such that 



« = K 



2tt 
d a N a 



m a — , V a. 
d 



(42) 



That is, the linear momentum of the emitted photon has 
to match the sum of three contributions: the momentum 
of the incident laser, kL; the initial momentum of the 
spin-wave, n2n/ (doN a ), and the contribution from the 
lattice periodicity, which enters through the reciprocal 
wavector m27r/do- In ID, condition (j4"2"|) has to be satis- 
fied only by the a — z component of these vectors, that 
is, only the projection of the momentum on the chain is 
conserved. In 3D, on the contrary, the equality has to 
be satisfied by all the vector components. The relation 
(|42p determines the maxima in the emission pattern de- 
pending on do and A, but it also determines the collective 
rates, through the normalization condition on I n (ty, see 
Eq. (00]). 

Since — 1 < Uq < 1, if A > 2do there is at most a 
single value of m, for which condition (|4"2"|) can be fulfilled 
for some vector uq. Thus, A > 2d is the directional 
regime in the emission of photons, whereas in the case 
A < 2do, we cannot ensure that photons are collimated 
in a single direction. In the following two subsections, we 



will study these regimes in ID and 3D. In particular, we 
will be interested in determining A9, the angular width 
of the photon beam in the directional regime, and r n , the 
collective emission rates. These quantities will be studied 
as a function of dimensionality, do, A, and N. 



B. Atom chains 

We consider for concreteness that kL points in the z 
direction, parallel to the chain axis. The vectors n, m 
are reduced now to scalars n, m, that correspond to the 
projections on the chain axis. Due to the symmetry of 
the problem, the photon distribution emitted by a spin- 
wave n depends only on the angle 9, through = cos 9, 



IJSI) 



A 



)N/2) 



1 f 1 sin 2 ({k L d (u z n - 1) 

47rT n N sin 2 ((fc L doK-l)-nf )/2) 

Taking the limit N ^> 1 and using 



/„(fi) = ^r/jr 1 ^) 

m 



(43) 



(44) 



k L d Q z . 7T 
— — (Ur) — 1) — n— 



The peaks in the photon distribution correspond to the 
emission of photons such that the z component of mo- 
mentum is conserved, 



A ii A 

Uo — 1 = m, 

" d N d 



(45) 



that is, there is a maximum whenever there is a value of 
to which satisfies this relation. In the directional regime, 
A > 2g?o, there is at most a single value of m which satis- 
fies (|45l) . However, in this case, momentum conservation 
only determines the value of 9 at the emission maximum, 
which implies that, in general, photons are emitted in 
cones spanned by different values of cj>. Only when the 
maximum happens at 9 — 0, or 9 — n, photons are colli- 
mated in the forward- or backward-scattering directions, 
respectively (see Fig. [2]). 

Let us study first the case 77 = 0, that is, the 
emission properties of the completely symmetric state. 
The forward-scattering contribution, I Q ^ , has an angular 
width given by A6>i£> = l/^k-^doN . The contribution to 
the emission pattern from each of the Bragg terms is 



Po 



Pa 



dni l ° ] (n) 



r a 



2 Po U \ 777 ^0. 



(46) 



Condition (|45|) leads to the result that the number of 
emission cones with 777 7^ is given by int(2d /A). A 
calculation made without resorting to the planewave ap- 
proximation, yields the same results, even for relatively 
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k" = K n ~ + k L ~ + m n/d 



k = K„ + k L 



FIG. 2: Photon emission by a spin- wave with linear momen- 
tum K n in an atomic chain. The maxima in the photon dis- 
tribution are at angles 9 such that momentum along the chain 
is conserved, up to reciprocal lattice wave vectors. 



small atom numbers (N — 20), see Fig. [31 However, 
the shape of the photon angular distribution in the ex- 
act calculation shows a departure from the sinc 2 -shape 
predicted by Eq. (|45|) . see Fig. [4j By using the normal- 
ization condition for /o(^)i we determine the probability 
of emission in the forward-scattering direction, 



1 



Pb = 



1 + 2 int 



'Mi 

v A 



(47) 



Thus, in the directional regime, A > 2do, all the diffrac- 
tion peaks but the forward-scattering one are suppressed. 
The collective emission rate is given by 



£0 

r 



A A . f2d 
H int — - 

Ad 2d V ^ 



(48) 
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1 

1 - cos(9) 



1.5 



FIG. 3: Angular photon distribution of the photonic mode 
emitted by the completly symmetric spin-wave in a chain of 
N = 20 atoms. We have used the exact egenvalues of the den- 
sity matrix, and followed the method described in the previous 
section. 

Finally, we calculate the emission rate for all the ID 
spin-wave states. They can be written as an integration 
over the contributions coming from different diffraction 




0.1 0.2 0.3 0.4 
1 - cos(0) 



FIG. 4: Angular photon distribution of the photonic mode 
emitted by the completly symmetric spin-wave from a chain 
of N = 20 atoms, A > 5do, calculated assuming periodic 
boundary conditions (continuous line), and performing and 
exact diagonalization of the master equation (dashed line). 
The distribution in the exact case shows a departure form 
the sinc 2 -shape due to collective effects. 



peaks, through the normalization condition, we have 

k L d 



r 



-E 

An ^ 



dxf_ 



■V 



-(x — 1) — n — 

v ' N 



The collective rates are thus determined by the number 
of diffraction peaks which appear in the emission pattern 
of each collective state n. Considering the limit in which 
k-^d^N ^> 1, that is, the delta limit for each of these 
peaks, we get the following distribution of emission rates, 

r» = 
r 

Xu> (W-^- + ^ + mt (- + ^j' (50) 

with xid — X/(2do), and 6{x) is the Heaviside function. 
This expression describes quite well the emission rates 
calculated without assuming periodic boundary condi- 
tions, that is, by diagonalizing the matrix Jjj given by 
([13)). see Fig. [5] In the limit d 3> A, we recover T n = f . 
On the contrary if A is comparable to do, emission rates 
are renormalized. In the directional regime, A > 2do, 
some of the states have an enhanced rate, T n — xid (su- 
perradiant), whereas there are states for which r„ = 
(subradiant) [l7|. Due to the relation between the pho- 
ton distribution and the emission rates l[4"9]). subradiant 
states correspond to spin-waves whose emission pattern 
is not peaked at a given value of 9. On the contrary, su- 
perradiant states are spin-waves whose emission pattern 
does contain a maximum as a function of 9. Note that 
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Eq. (|50|) does not describe the case of Dicke superra- 
diance, which would predict a single superradiant state, 
with To = AT. The reason is that we have assumed 
condition Nk^do 1, that is, where are always in the 
regime in which the light wavelength is much smaller than 
the size of the chain. Our results for ID are consistent 
with Ref. [l8(, where the collective light emission form 
an atomic chain was studied with the quantum jump for- 
malism. 




X = 4 d 



o 
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FIG. 5: Calculation of the collective rates of a chain with N = 
20 atoms. Filled squares: result predicted by the planewave 
approximation (|50p . Empty circles: eigenvalues of the matrix 
Jij defined by (|14|l . In the last case, to each eigenvalue of Jij 
we assign a wavenumber n that corresponds to the maximum 
Fourier component of the corresponding eigenvector. 



To summarize the situation in ID, in the directional 
regime A > 2do, superradiant spin- waves emit photons 
into a single emission cone, with a rate that is enhanced 
by Xid- This effect can be used to generate photons that 
are collimated along the chain axis. Note that, in general, 
the atom-photon mapping induced with the A-scheme of 
Fig. [T]may compete with other radiative processes, such 
as the radiative decay from |a) back to |s), or the radia- 
tive decay from \a) to other atomic levels that are not 
included in the A-scheme. This problem can be solved 
by enhancing the atom-photon mapping rate, choosing 
Xid 1. Since \id does not depend on the chain size, 
this implies to choose A ^> do . Another way out of this 
problem is to use a cycling transition (see Appendix [D| . 



C. 3D Atom lattices 

Contrary to the ID case, in 3D it is not simple to ob- 
tain closed expressions to describe the emission in the 
planewave approximation. In order to get a simpler pic- 
ture, we replace the sinc-function in the definition 
by a gaussian which is normalized in the same way, 



f N (x) w s/ttN x 



This approximation is justified in the limit in which the 
diffraction peaks are narrow enough, such that they do 
not overlap, that is, N a khdo ^> 1. Assuming that the 
length of the lattice is the same in any spatial direction, 
N x = N y = N z , then the angular photon distribution for 



a state n is given by 

j n (n) = E J n ml ( fi )' 

m 

n 1 ' 4vrr n 



(52) 



A': 



A given term 7n m '(f2), has a non- negligible contribution 
only if 



|k L do + 2nn/N + 27rm| 
khdo 



1. 



(53) 



since otherwise, there are no photons which satisfy the 
energy-momentum conservation. 

Let us study first the directional regime, A > 2do- In 
this case we find the following two possible situations: 

(i) Superradiant states - Spin-waves for which there 
exists a value of m, say m c , such that ([53"]) is satisfied. 
The condition A > 2do ensures that there is a single value 
m c , such that the photon distribution can be simplified, 



Z n (0) 



i(n). 



(54) 



Thus, there is a single emission peak in the direction, 
k L d + 2irn/N + 27rm c 



(55) 



Photons are collimated in a beam with width A03D = 
l/(khdoN 1 ^ 3 ). From (|54|) . and the normalization condi- 
tion, we can determine the emission rate, 



r 



N 

47T 



two K L"0 L 

dtt 7T 3 / 2 e I 2 nn 2 



J^,„. ilL_!il — Tr-j^-— 7rm c ^ 



Jn \ 2do 



X3D- 



(56) 



(51) 



Thus, superradiant states have a decay rate that is en- 
hanced by the optical thickness, X3D- 

(ii) Subradiant states - Spin-waves for which there is 
no value of m which satisfies ([5"3"|) . In this case, there 
is not a single dominant contribution m, such that, di- 
rectionality in the emission of photons is not guaranteed. 
Also, the normalization condition on leads directly 
to the result r n = (T7j . 

From this analysis, we conclude that superradiant spin- 
waves are interesting in the context of quantum infor- 
mation, since they emit collimated photons with an en- 
hanced emission rate. However the emission direction 
is determined by a value of m c which has to be cal- 
culated for each particular spin-wave. There are two 
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relevant cases in which this situation becomes simpler, 
since energy-momentum conservation ensures that the 
only Bragg contribution is m c = 0. The first one is 
the case A > 4do, because 



kg I 



< 2kr 



.n a 2ir 2-7T 1 2n ... 

l NT +m ^ > 2T Q (lf m ^ 0) 



(57) 



and thus for values A > 4d the energy-momentum con- 
servation condition (|4"2"f is only fulfilled by m c = 0. Sec- 
ond, assume that the spin- wave is created by the absorp- 
tion of a photon within the A-configuration of Fig. Q] 
such that the spin-wave linear momentum can be writ- 
ten like K n = q — kL, with q = k^. Then, the energy- 
momentum condition is now 



q + k = m 



d ' 



(58) 



which is satisfied only by m c = if A > 2d,Q. Indeed, 
the spin-wave created in this way is always superradi- 
ant, since K n automatically satisfies Eq. (|53p. A similar 
result is obtained in Ref. [l9( for the case of a cloud 
of atoms, where the rate of spontaneous emission of a 
state created by the absorption of a photon, is shown to 
increase with the optical thickness. 

To summarize the situation in 3D, in the directional 
regime A > 2do, superradiant spin- waves emit collimated 
photons with a rate enhanced by X3D- Contrary to the 
situation found in ID, in 3D lattices X3D oc N x , and thus 
it can be increased by increasing the size of the system. 
In this way, the atom-photon mapping could be much 
faster than other competing decay channels. 



IV. EMISSION OF LIGHT BY HOT ATOMIC 
ENSEMBLES 

The formalism presented in the previous section can 
be extended to situations in which atoms are not at fixed 
positions in space. In this case, one expects the photonic 
state not to be a pure state, but a mixed one. This 
is relevant to the description of hot atomic vapors. We 
follow the discussion presented in subsection III D[ and 
consider the case of fast motion (r mot <C I\A The results 
presented here are in agreement with Ref. [201 ] , where the 
photonic mode emitted by an atomic ensemble is studied. 
Our first task is to reformulate the problem of photon 
emission by finding the master equation that describes 
this situation. 



A. Master Equation for atomic ensembles 

In the fast motion limit, particle positions may be de- 
scribed as a set of independent random variables, {ij}. 



We choose for simplicity a gaussian distribution proba- 
bility, 



1 



n 3/2 L 3 



-(r/Lf 



(59) 



The particle positions fulfill the following identity, which 
will turn out to be the basis for the following calculations, 



' ) Uot = 6 jl + (l-6 jl )e-« 2L2 / 2 , 



(60) 



where (. . . ) mo t is an average over the atomic positions. 

To calculate the coefficients of the master equation we 
start from Eq. (|12p and perform the following average, 



J i ■ 



k 



g k 5(uj k -uj h ){e 



(k-k L )( ri -r 3 



(61) 



were we have neglected the Cauchy principal value con- 
tribution, since it leads to an energy shift that does not 
play any role in the discussion that follows. After using 
(|6H)) . integrating over k, and taking the limit k^L ^> 1, 
we get 



Ju = r/2, 



Ji 



4(fc L L) 



(62) 



One can readly diagonalize the matrix and obtain 
the eigenspaces of the master equation. The first is the 
completely symmetric state (n = 0), 



M j0 = 



1 



N 



Jo 



(Xc 



1) 



(63) 



where Xcn — (N - l)/(2(fc L L) 2 ) is the optical thickness 
of the atomic ensemble. Note that in the limit N ^> 1, 
and defining the atom density, n at = N/L 3 , the optical 
thickness can be recast in the more familiar form Xcn = 
n a t X 2 L/2. The second eigenspace is spanned by the spin- 
waves orthogonal to Aijo, 



3 



Mjn =0, J n = 



1 - 



1 



2(k L Ly 



if n ^ 0. 



(64) 

Collective effects happen if Xcn ^> 1 . In this case, there is 
a single superradiant spin-wave mode, corresponding to 
the completely symmetric state, and N — 1 states which 
decay with the single atom emission rate, J n = r/2. 



B. Single photon state 

The spin-wave mode n = is the only one to show 
collective effects, and it is the collective state which can 
be created in experiments with atomic vapors. For these 
reasons we focus on the following on its properties. First, 
we consider the emission of light by the initial atomic 
state with a single excitation, 



l*o> 



1 



N ^ 3 



(65) 
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The atom-photon mapping can be extended to this situa- 
tion by considering first that the atomic state is mapped 
into a given photonic state, and then by performing the 
average on the atomic positions. This gives as a result a 
photon density matrix, 

l*0)at|0) ph (0| -> |0W, 

P = ^Pkk'akl°)ph(0|ak'- (66) 

k,k' 

Note that the mapping is now form pure to mixed 
states. Unfortunately, with the statistical properties of 
the atomic positions considered here, it is not possible to 
ensure that this is still a gaussian map. We can, however, 
use Eqs. (|24I60|) to get 



Pkk' 



30,k<Po,k'/mot 
1 



(1 ~ e)0o,k$j, k ' + e P> 



1+Xc 



(67) 



We work, for clarity, in the interaction picture with re- 
spect to -ffi m . The coherent component of the photon 
density matrix is the pure state |0o)ph, 



1 



'",k = \/H y/N- lgk— f A , T , 



(68) 



whereas p is a normalized mixed state which describes 
the incoherent (isotropic) component, 



Pkk' 



(1 +Xen)9k9k'e 



-|k-k'| 2 L 2 /4 



(-i(u>k - wl) + Jo) (i(uk' - ^l) + Jo) ' 



(69) 



Under the condition Xcn S> 1, the main contribution to 
the photon density matrix is the pure one, |0o)ph- The 
parameters which describe the directionality in the pho- 
ton emission can be readly evaluated. The emission pat- 
tern has both coherent and incoherent contributions, 

/o(fi) = r oh (n) + i^(n), 
r oh (n) = ( i_ e )^l e -(fe L u fi -k L) 2 L 2 /2 ; 

47TXcn 



2T(n) = e-. 



(70) 



Thus, photons are collimated in the forward-scattering 
direction, with angular width A8 cn = l/(ki J L). The devi- 
ation of the photon emission from directionality is deter- 
mined by the probability of emission out of the forward- 
scattering cone, £ = J G?f2/ nc (f2) = e. 

Finally, we calculate the purity of the photon state 
which is defined as 

V = Tr(p 2 ) 

= (1 - e) 2 + e 2 Tr (p 2 ) + 2e(l - e)^o|/#o>. (71) 

This quantity is very relevant when using photons in 
quantum information processing, since it describes the 
efficiency in the process of interference of photonic modes 



in a beam-splitter, being V = 1, the case corresponding 
to a pure state. To calculate V, we notice first that 



k,k' 



2(fc L L) 



(72) 



Also because of the Cauchy-Schwartz inequality, this im- 
plies that 

|(0o|p|^)| 2 <Tr(p 2 ), (73) 

such that, 

V = Tr(p 2 ) = (1 - e) 2 + 0(1/*lL) « 1 - 2/ Xcn . (74) 

Thus, in the limit k^L 3> 1, the purity of the photon 
state is solely determined by X3-D- 

C. Multiphoton state 

We study now the properties of the multiphoton case. 
The mapping is from a state with M atomic excitations in 
the completely symmetric state to a multiphoton mixed 
state, 



|0>at|0)ph-|<W (75) 



We assume the low-excitation limit (M -C iV). In Ap- 
pendix [B] we show that the mixed photonic state can be 
written like 




p = (l-e) M |$) ph ($|+p, 

M 



7=fe^ ka k) l°>ph- (76) 



That is, it is a sum of a pure state consisting of M pho- 
tons in the photonic mode defined by ([55)1 . and the (not 
normalized) mixed state p. In order to determine the 
purity of p ph we should study the mixed contribution, 
which is far more complicated than in the single photon 
case. Insted, we notice that 



V = Tr(p 2 ) > (1-e) 



2M 



(77) 



which allows us to obtain a lower bound for the purity. 



V. IMPLEMENTATIONS 

A. Collective light emission in ultracold atoms in 
optical lattices 

Ultracold atoms in optical lattices are an ideal sys- 
tem for the observation of the effects described in this 
work. For example, atoms in a Mott phase [2l[ would 
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be ideally suited to study collective light emission from 
a square lattice. In this setup, atoms are placed at dis- 
tances that are comparable to optical wavelengths, since 
potential wells in a standing-wave are indeed separated 
by do = A sw /2, with A sw , the wavelength of the coun- 
tcrpropagating lasers that create the lattice. The con- 
ditions for the directional regime are met by using an 
optical transition such that A > A sw . Under this con- 
ditions, ultracold atoms in optical lattices are ideal to 
form an atom- light quantum interface [Til . I22I |23| , where 
the quantum state of light can be indeed manipulated, 
as show recently in [2J]. The properties of the emitted 
light may also be used to measure properties of quantum 
many-body phases in optical lattices [25| . Our theory can 
also be applied to the collective matter-wave emission in 
optical lattices proposed in Ref. 26]. 

We propose now an experiment in which the renor- 
malization of the collective rates may be observed by 
performing the following steps (see Fig. [6]): 



(i) 



(") 




\g) \s) 



„/ \ / \ 

' K = k LA -k LB k = K +k L K = R l 



A L,B 



FIG. 6: (i) Creation of a spin- wave with linear momentum 
\iL.A + ^l,b in an optical lattice, (ii) Mapping from the spin- 
wave to a photonic mode. 

(i) First, we use a protocol to initialize the atomic 
quantum state of the lattice. Assume that all the atoms 
are initially in the ground state \g). The A-scheme of Fig. 
[I] may be used for the initialization by shining the lattice 
during a short time T with two lasers with wavevectors 
ki^Aj kL,B, and Rabi frequencies J7l,a, ^l,b, and de- 
tuning A, such that Ol.a, ^l,b <C A. Under these con- 
ditions, a coherent state of spin-waves with momentum 
K = kL,A ~ kL b is created, 



|*) at ex e-^x/2)6+ |0)s 

b * = -^E e * Krj ^> 



(78) 



where ^ab = ^la^l,b/A, and condition QabT <C 1 
has to hold in order to ensure that the system is in the 
low excitation regime. 



(ii) After the creation of the collective atomic state, 
this can be realeased by a second laser with momentum 
kL in the A-scheme considered along this work. The anal- 
ysis presented in section UTT1 can be applied to study the 
emission of photons by setting K n = kL,A — kx,B- For 
example, if kL = k^B, then the spin- wave emits a photon 
in the direction kL,Ai with a superradiant rate X3£>r. On 
the contrary, if K n does not satisfy the condition (|53p . 
then the emission rate will be suppressed, up to finite 
size effects. 



Deterministic generation of entangled states of 
photons 



An application of our ideas to generate photons in a 
deterministic way requires a system where experimen- 
talists are both able to reach the directional regime of 
collective light emission, and to initialize the atomic sys- 
tem in a given spin-wave state. This idea may find use- 
ful applications in quantum cryptography [2^1. quantum 
computation [28j |. and quantum litographv |29l|. 

Before going into the description of particular experi- 
mental setups, we discuss how to extend our formalism 
to generate photons that are entangled in polarization. 
For this, we consider the double- A scheme of Fig. 
in which there are two excited levels, |si), \s2), and a 
ground state level \g). The atom-photon mapping is de- 
scribed in the same way, by defining the corresponding 



Holstein-Primakoff operators b 



Assume that the 



l s i), l s 2) states decay by emitting photons with differ- 
ent polarizations and creation operators a]*~ k , aj k - If the 
conditions for momentum conservation explained along 
this work are fulfilled, then entangled spin waves will be 
mapped into entangled photons in polarization. For ex- 
ample, consider the following initial spin-wave state, 



l*)at = 



— ( b + b + + b + b + 



(79) 



where we have used the notation b 



(VViv)E,- 



JKr, 



b+. 



After the atom-photon map- 



ping, the following photonic state is created, 



|0>) 



P h 



V2 



a 



+ 



l,K B +k L u 2,K y i+k L 



) P h.(80) 



In this way, the ability to generate entangled spin-waves 
is equivalent, by virtue of the atom-photon mapping, to 
the ability to generate entangled photonic states. We 
propose three experimental set-ups, where this is possi- 
ble and discuss the conditions that are required for the 
implementation of this idea. 
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FIG. 7: Double A-scheme for the generation of photons en- 
tangled in polarization. 



1. Trapped ions 

This system is specially appealing from the point of 
view of the creation of entangled states, since quantum 
gates allow us to create deterministically any collective 
state. This idea has been demonstrated in the creation of 
the K = (W) state, in Q. The first issue that we have 
to deal with, is whether conditions for forward-scattering 
of photons hold here. Usually ions are arranged in chains, 
with the peculiarity that the distance between ions is 
not constant. However, if condition <jp v < A/2, with <ig v 
the average distance, we still get light emission in the 
forward-scattering cone only. We show this fact by per- 
forming a calculation of the emission pattern of a single 
atomic excitation in the K = state, using the the- 
ory presented in the previous sections, see Fig. [5) Note 
that directionality is achieved even for a relatively small 
number of ions (N = 10). The main difficulty for the im- 
plementation of this idea with ions lies on the fact that 
ion-ion distances are usually in the range of a few /xm, 
and thus condition dg v < A/2 is not fulfilled when consid- 
ering optical wavelengths [301 ■ A way out of this prob- 
lem is to use optical transitions which lie in the range 
of A > 5^m, like for example A( 2 D 3 / 2 - 2 Pi/2) = 10.8 
^m in Hg+, or A( 2 D 3/2 - 2 D 5/2 ) = 12.5 ^m in Ba+. A 
particularly interesting initial atomic state is the spin- 
wave ([75)1 with = 0, and k# = — 2kL, (all vectors in 
the direction of the atomic chain), since it leads to the 
emission of an entangled photon pair in the forward and 
backward-emission directions. 




1 - cos(0) 

FIG. 8: Angular photon distribution of the photonic mode 
emitted by the completely symmetric spin-wave in an ion 
chain with N = 10 ions. 



More efficiently, one could avoid the use of quantum 
gates by using the dipole-blockade mechanism with Ry- 
dberg atoms, which allows us to generate W-states, as 
well as states which emit Fock states with a number M of 
photons 0. Interactions between excited atomic states, 
like those that take place in Rydberg atoms, can be also 
used to generate photons entangled in polarization. This 
can be achieved in a single experimental step, without 
the need for quantum gates, if the proper configuration 
of atomic interactions is chosen. As an example, consider 
the level configuration shown in Fig. [51 and interactions 
between excited states such that atoms in levels |si), |s 2 ), 
interact strongly only if they are in the same excited 
state, that is, U\\ = t/ 22 = U, but £/ 12 = 0. We ap- 
ply two lasers with wavectors kA,B and Rabi frequencies 
^a.b, detuned with respect to the \g) - \s a .b) transition, 
such that Aa = — Ab = A. If condition Aa.b 3> ^a,b 
is fulfilled, then the lasers induce a two-photon transi- 
tion with Rabi frequency f2 e gp = JIa^b/A. Furthermore, 
if f2 c ff -C U, states with two atoms in the same excited 
state are not populated. Under these conditions there 
are two possible excitation channels, depicted in Fig. [SI 
which give rise to the linear combination (|79[) . 



2. Ultracold atoms in optical lattices 

We have shown above that optical lattices are well 
suited to reach the regime of directionallity in the emis- 
sion of photons. The main issue here, contrary to the case 
of trapped ions, is to find a way to create efficiently the 
initial spin-wave states in a deterministic way. Although 
one could think of peforming quantum gates between ul- 
tracold neutral atoms to generate collective atomic states 
[3lL [3^ | , this procedure faces the difficulties of quantum 
computation in this system, like for example, how to 
achieve single atom addressability. 



3. Atomic ensembles at room temperature 

The very same techniques which can be applied to Ry- 
dberg atoms in an optical lattice can also be used in the 
case of hot ensembles. On the one hand, this setup has 
the advantage that atoms do not need to be cooled and 
placed in an optical lattice. On the other hand, it can 
be described by a statistical distribution of particles, and 
thus suffers from the fact that high efficiency in the re- 
lease of photons is achieved under more severe conditions 
of particle density and atom number, as discussed above. 
However, densities which are high enough to fulfill the 
requirement £ <C 1 have been recently reported in [33| . 
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\g)ig)j 

FIG. 9: Scheme for the generation of entangled spin-waves by 
using the Rydberg blockade. 



C. Photon up-conversion 



Finally, the ideas presented in this work find an intere- 
seting application in the efficient up-conversion of pho- 
tons. For example, transitions between the states \g) and 
|s), are typically in the range of GHz. Application of the 
A scheme of Fig. [T] provides a way for up-conversion of 
the microwave photon to an optical photon. Also, by us- 
ing several lasers, one could induce TV-photon transitions 
from \g) to |a), and emit a single high frequency photon. 
By using the ideas presented along this work, this could 
be done in such a way that the up-converted photons are 
collimated and, thus, they can be efficiently collected. 



VI. ACKNOWLEDGMENTS 



We acknowledge interesting discussions with M. 
Lewenstein. Work supported by the E.U. (SCALA), and 
the D.F.G. through the excellence cluster Munich Ad- 
vance Photonics (MAP). 



APPENDIX A: VALIDITY OF NEGLECTING 
PROPAGATION EFFECTS 



When performing the integration in k in Eq. (j2"Tj) , we 
have to deal with the following factor inside the integra- 
tion, 



e ,i(k-k L )(r j -r J ,) 



(-i{uj k - Wl) - Jn) - w l) - Jn') ' 



(Al) 



which can be evaluated in the following way, 
1 

(-i(aj k - w L ) - J*) (i{u k - uj l ) - J n >) 
1 1 



-i((Jk - Wl) - J'n i{Uk - W L ) - Jn' 

w -27T(5(wfc - uj L ) * 
where he have aproximated 



Re 



1 



-i(uk - w L - A„) - T„/2 



ir5{uj k -w L ).(A2) 



This approximation is justified by the following argu- 
ment. The left-hand side of (|A2|) is a function of k with 
width r„/c. The latter has to be compared with the 
width of the exponential in (|A1[) . which is, roughly 1/L. 
Thus, the approximation holds in the limit T n /c -C 1/L, 
that is, whenever the emission time is shorter than the 
propagation of the photon through the sample, some- 
thing that is well justified in the trapping setups consid- 
ered in this work. The energy shift can also be safely 
neglected, since condition A„ < wj, is well justified by 
the atomic transitions used in cold atom setups. 



APPENDIX B: MULTIPHOTON STATE 
EMITTED BY AN ATOMIC ENSEMBLE 

Our starting point is the atom-photon mapping defined 
in Eq. (|75[) . By applying our method we get 



1 1 

Ml 'N 1 ' 1 ' 



E 

3M 



e i(k-k L )rj 



i(u> k - ujl) - r /2 



|0) P h(0| ( 



i(k'-k L )r • 



-i(wk> - ujl) - r /2 




(Bl) 



After expanding the two parenthesis to the power of M 
we can reexpress this equation as a sum of 2M products 
of exponentials of the argument z(k — kjjrj. There are 
Gm) terms such that all the random variables are differ- 
ent. These terms can be summed up to a contribution 
which yields a pure state, such that p is the sum of a 
pure and a mixed state, 



P = Ppurc + P, 



(B2) 



where p pur o, P, are not normalized. The pure contribu- 
tion has the form, 



Ppu 



|0) ph (0| Yl 




M 



9ke 



^ i(u k - uj l ) - r /2 

\ M 



i(uj k , - Lu h ) - r /2 



Ok' 



(B3) 
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which is a pure state with M photons. To normalize this 
state, we use the limit (J^) /N M = N AI , (N > M), and 
obtain the final result (f7fj|) . 



APPENDIX C: 3D LATTICE, LIMIT A < do 

In the case of atoms in a 3D square lattice, in the limit 
A <C do, there are in principle many diffraction peaks 
which contribute to the emission pattern. In the fol- 
lowing we show, that, at least for the completely sym- 
metric state, if the size of the system is large enough, 
the forward-scattering contribution is the most impor- 
tant one. 

Recall the definition of the probability that the spin- 
wave emits a photon in the m Bragg scattering peak, 
Po™' PI]) - The forward-scattering contribution m = 0, is 
given by 



[o] r 

Po = T^X3D- 
1 



(CI) 



We consider for concreteness Icl = k^z. To determine the 
probability of scattering out of the forward-direction, £ 
we calculate the sum over the contributions with m^O, 

£ = E^ m] = 

m/0 



L V — [dSl ^3/2 e -(« fi ---^m) 2 (fe L d ^/2) i 

To ^ n 4tt ' 



• (C2) 



Considering the limit ki^doNx/2 ^> 1, the angular inte- 
gral yields, 

£ = - ^t 3/2 — c -(i-\i r ™+m 2 (k Ij d N !C /2) 2 



r (k L d N x y 



(C3) 

In the limit do 3> A, we can replace the sum by an inte- 
gration over m, and get the result, 

£ = i-> (C4) 

l o 

Together with the normalization condition ^2 m p m = T 
this leads to 

£ = —^—, ^- = l + X3D- (C5) 

i + X3D r 

Thus, as long as the optical thickness is large, the com- 
pletely symmetric state is superradiant, and photons are 
collimated in the forward-scattering direction. Note that 
the regime A <C e?o is very different than the directional 



regime, in the sense that there is always a non-zero prob- 
ability of emission out of the forward-scattering cone. In 
the fully directional regime (A > 2do), on the contrary 
£ = 0, up to finite size effects. 

APPENDIX D: CYCLING TRANSITION 

Up to now we have neglected the decay form the aux- 
iliary level | a) back to the excited state |s). This process 
can be suppressed by considering, for example, a level 
configuration like the one presented in Fig. 1101 |s) and 




FIG. 10: Level scheme for the implementation of the atom- 
photon mapping, avoiding incoherent processes. 

\g) could be two hyperfine states of the electronic ground 
state manifold, and \a) be chosen such that the |a)-|g) 
is a cycling transition. For example, these levels could 
be found in 87 Rb: \s) -> |S 1/2 ,F = l,m F = 1), \g) -> 
|S 1/2 ,F = 2,m F = 2), and \a) -> \P 1/2 ,F = 3,m F = 3). 
Spontaneous emission from \a) to |s) is forbidden by se- 
lection rules on m F , but atoms decay from \a) to \g) 
with a rate r ga . Fields with Rabi frequencies fii, f2 2 , 
and detunings Ai, A 2 , are coupled to the transitions \s)- 
\g), \g)-\s), respectively. Under some conditions, lasers 
induce a two-photon transition that is equivalent to a 
A-configuration such as the one considered in this work 
with = ^eff = Oif^/Ai. 
We consider two situations: 

(i) A 2 3> r ga . The atom-photon mapping is possible 
under conditions: fl e g <C A 2 (adiabatic elimination of 
the level |a)), y/Nfli <C Ai (to avoid real transitions 
from |s) to \g)), and r ga (f2 cff /A 2 ) 2 » NT ga (Q 2 /A) 2 , 
or equivalently, £!i/A 2 ^> >/~N, to avoid the dephasing 
induced by the transfer of atoms from \g) to \a). 

(ii) A 2 = (resonant case). Conditions: fi c ff <C r ga 
(adiabatic elimination of \a)). Q1/A2 3> V^V, and 
r ga (fi c ff/r ga ) 2 > iVr ga (f7 2 /Ai), or equivalent^, flj > 
iVr ga (same reasons as (i)). 
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